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I.  Introduction 


The  M4S3A1  projectile  was  the  result  of  a  modification  of  the  XM4S3  configuration. 
This  modification  consisted  of  shortening  the  boattail  length  from  0.5  calibers  to  0.25 
calibers.  This  was  determined  to  be  neccessary  as  the  XM4S3  exhibited  a  very  large 
positive  Magnus  moment  at  Mach  numbers  less  than  1.0.  This  large  positive  Magnus 
moment  caused  the  projectile  to  be  dynamically  unstable  (fast  mode)  at  fairly  low  angles 
of  attack  (yaw  angles)  resulting  in  projectiles  falling  short  of  their  design  range.  Reference 
1  presents  the  analysis  of  spark-range  and  yawsonde  firings  of  the  XM4S3.  Also  presented 
in  Reference  1  are  the  results  of  yawsonde  firings  of  a  modified  XM4S3  (Mod  21B)  which 
is  now  the  configuration  of  the  M4S3A1. 

Figure  1  is  a  sketch  of  the  M4S3A1.  The  M483A1  (20  each)  was  fired  in  the  Ballistic 
Research  Laboratory  (BRL)  Transonic  Range  (TR)  Facility  in  the  mid  1970’s  along  with 
105mm  and  155mm  models  of  the  XM795.  The  results  of  these  tests  are  given  in  Reference 
2.  Non-linear  Magnus  moment  data  are  obtained  from  reductions  utilizing  the  Eglin  Six- 
Degree-of-Freedom  (6DOF)  Aeroballistic  Range  Facility  Data  Analysis  System  (ARFDAS) 
program.  Reference  2  discusses  the  very  non-linear  Magnus  moment  discovered  during 
initial  attempts  utilizing  cubic  and  quintic  expansion  (with  angle  of  attack)  terms.  The 
6DOF  model  was  changed  to  perform  a  table  lookup  of  estimated  values  of  the  Magnus 
moment  derivative  ( Cnpa ).  This  change  produced  satisfactory  analysis  results.  However, 
with  only  20  projectiles  tested  between  Mach  numbers  0.65  and  1.91,  the  many  data  gaps 
in  Mach  number  and  angle  of  attack  resulted  in  a  partial  aerodynamic  data  package. 

Recently,  the  exterior  shape  of  the  M4S3A1  has  been  selected  as  a  baseline  for  the 
development  of  several  new  projectiles.  The  aerodynamic  data  base  of  the  M483A1  was 
determined  not  to  be  of  sufficient  quality  when  the  importance  of  development  programs 
is  considered.  The  BRL  designed  a  TR  test  program,  consisting  of  45  additional  firings, 
to  supplement  the  1970’s  data.  This  program  was  carried  out  from  1987  thru  1989.  This 
report  discusses  the  results  of  these  new  tests  in  combination  with  tne  earlier  tests  which 
were  re-analyzed  in  concert  with  the  new  test  data. 


II.  Procedure 

1.  Test  Facility 

The  projectiles  were  fired  thru  the  207  meter  long,  spark-shadowgraph  instrumented 
section  of  TR.  Reference  3  describes  the  details  of  this  facility.  The  shadowgraphs  taken 
are  then  used  to  determine  the  position  and  attitude  of  the  test  projectile  as  a  function  of 
time  and  distance.  There  are  twenty-five  shadowgraph  stations  in  the  range.  During  the 
tests  conducted  in  1975,  only  eighteen  of  the  twenty  five  stations  were  timed. 

The  expected  measurement  accuracy  is  0.1  degrees  in  attitude  and  0.003  meters  in 
position. 
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Table  1.  M483A1  Physical  Properties 


Period 

Diameter 

Leugth 

CG 

CG 

Mass 

mam 

H 

mm 

mm 

mm 

Cal. 

kg 

E9B31 

1975 

154.74 

3.65 

■jityJI 

1.695 

1988 

154.74 

896.7 

|  562.9  1 

3.64 

iHil 

wmm 

1.687 

2.  Test  Projectiles 

The  projectile  configuration  is  shown  in  Figure  1.  Actual  projectile  hardware  was 
utilized  for  the  tests.  The  dimensional  characteiistics  and  physical  properties  of  each 
projectile  tested  were  measured  prior  to  firing.  Key  characteristics  are  given  in  Table  1. 

As  shown  in  Table  1,  only  minor  differences  between  the  1975  and  1988  projectiles 
were  measured.  The  reference  CG  for  the  data  presented  in  this  report  is  3.64  calibers 
measured  from  the  nose.  The  stability  computations  utilize  the  1988  physical  properties. 


3.  Data  Reduction 

Data  reduction  procedures  were  formulated  by  using  epicycles  to  fit  the  spark-range 
data  (Ref  4).  This  method  is  often  called  the  linear  theory.  The  data  analysis  technique 
and  system  utilized  for  this  report  is  fully  described  in  Reference  5.  Slight  modifications  to 
the  baseline  system  (Eglin-ARFDAS)  were  made  for  compatibility  with  in  place  procedures 
at  BRL.  Some  notational  and  scaling  differences  exist  between  References  4  and  5.  The 
system  used  in  this  report  allows  for  the  simultaneous  reduction  of  up  to  five  experimental 
data  sets.  This  is  a  powerfull  technique  as  the  only  additional  unknowns,  required  to  be 
determined,  are  the  initial  conditions  of  each  added  data  set.  The  aerodynamic  coefficients 
to  be  determined  are  common  to  each  data  set.  A  toted  of  sixty-five  single  shots  and 
nineteen  multiple  fit  groups  were  reduced  during  this  activity  utilizing  the  six-degree-of- 
freedom  methodology. 

The  analysis  process  is  judged  to  be  adequate  when  the  resultant  fits  to  the  data  set 
approach  the  measurement  capabilities  of  the  facility  presented  above. 

The  standard  method  of  modeling  the  expected  non-linearities  with  angle  of  attack 
for  each  coefficient  and  its  derivative  is  given  below: 

e  =  sin(a) 

a  =  Total  Angle  of  Attack 
Axial  Force  Coefficient 

Cx  =  CXo  +  CX2e2  +  Cx4e4 
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Normal  Force  Coefficient 


C/V  —  C^aot  +  C/Vo3f3  +  Ctfc  5C5 


Csa  ~  CNao  +  CnoZ*2  +  C^aS^ 


Magnus  Force  Coefficient 
C/Yp  =  Cypaot  "t"  CVpor3e 

Cypa  =  Cypao  ”1*  CVpo,3C 


Pitching  Moment  Coefficient 


Cm  —  Cmao£  “I"  Cma3 C  "t"  Cma s€ 


Cma  —  Cmao  1"  Cma3 £  "h  Cma$€ 


Pitch  Damping  Coefficient 


Cmq  —  Cmqo  +  Cmqa2e2  +  Cmqa  46* 


Magnus  Moment  Coefficient 


Cnp 


c 


npa 


Cnpao  *■  +  Cnpa3e3  +  CnpaS£ 


Cnpao  *4“  Cnpa3€  4"  Cnpa&t 
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A  lesson  learned  from  Reference  2  was  that  adding  additional  polynomial  terms  to 
account  for  M483A1  non-linearities  in  Magnus  moment  was  not  sufficient.  The  1975  tests 
contained  data  with  projectile  yaw  exceeding  fifteen  degrees.  The  technique  chosen  for 
the  analysis  contained  in  this  report  (different  from  Reference  2)  was  to  create  a  table 
of  Magnus  moment  (Cnp)  as  a  function  of  angle  of  attack  (a)  and  Mach  number.  As  the 
projectile  motion  is  numerically  integrated  in  the  6D0F  code,  a  two  way  table  interpolation 
is  done  to  determine  Cnp  at  each  angle  of  attack  and  Mach  number.  This  value  of  Cnp  is 
then  added  to  the  standard  Cnpa  equation  by  the  following  technique. 

Cnpo,batt  =  Cnp/e 
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now: 


Cnpa  —  Cnpab OJ€  4*  Cnpabalt  *  Cg  -f"  Cnpao  4"  Cr»por3^  4"  CnpaS^ 
where:  Cg  is  to  be  solved  for. 


ill.  Results  and  Discussion 


Each  test  projectile  was  initially  subjected  to  single  fit  reductions.  The  aerodynamic 
coefficients  of  the  projectiles  are  determined  by  individual  reductions.  Following  the  sin¬ 
gle  fits,  projectiles  with  similiar  Mach  numbers  are  grouped  together  for  a  multiple  fit 
reduction.  These  groups  usually  contain  the  extremes  of  mean  squared  yaw  (low  and  high) 
available  at  that  groups  Mach  number.  A  maximum  of  five  projectiles  can  currently  be 
utilized.  In  some  cases,  more  projectiles  (greater  than  five)  are  desired.  However,  this 
analyst  is  content  with  five  data  sets,  reducing  more  data  sets  is  mind  boggling.  The 
tabulated  single  fit  results  are  given  in  Table  2  and  the  multiple  fit  results  are  presented 
in  Table  3. 

The  precision  of  the  single  fits  are  seen  to  better  the  stated  range  accuracy.  This 
analyst  believes  that  the  accuracy  of  the  calibration  and  film  reading  has  been  considerably 
improved  in  recent  years.  The  quality  of  the  multiple  fit  validates  this  premise  in  that  they 
also  better  the  stated  accuracy. 

The  1975  rounds  can  be  identified  as  having  ID’s  of  T13XXX  and  T14XXX  while 
the  1988  series  have  ID’s  of  T30XXX.  The  multiple  fit  groups  contained  combinations  of 
rounds  fired  in  1975  and  1988.  The  rounds  contained  in  each  multiple  fit  are  designated 
in  the  left  hand  column  of  Table  3. 

The  major  coefficients  will  be  discussed  individually.  Emphasis  has  been  placed  on 
the  multiple  fit  results.  Only  minimal  effort  was  expended  in  attempts  to  determine  non- 
linearities  from  single  fits. 


1.  Axial  Force  Coefficient 

The  zero-yaw  axial  force  coefficient  (also  Cd0 )  is  plotted  versus  Mach  number  in  Fig¬ 
ure  2.  The  single  and  mutiple  fit  are  separately  graphed.  Little  difference,  if  any,  is  noted 
when  comparisons  to  Reference  2  are  made.  The  error  in  the  axial  force  coefficient  at  zero 
yaw  is  estimated  to  be  less  than  0.004.  The  yaw-axial  force  (Cxi)  is  shown  for  each  mul¬ 
tiple  fit.  When  the  average  yaw  levels  were  under  5  degrees,  this  coefficient  was  estimated 
and  held  constant  during  the  analysis.  Reference  2  provided  estimated  values  for  C.\'2  as 
the  average  yaw  levels  were  large  for  the  M483A1  and  similiar  XM795. 
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2.  Normal  Force  Coefficient 


The  normal  force  coefficient  determination  is  dependant  on  the  magnitude  of  projectile 
yaw.  Assuming  a  typical  gun  launch,  with  a  first  maximum  yaw  of  3  degrees,  the  size  of 
the  slow  arm  A's  will  be  1.5  degrees  For  the  M483A1,  this  will  result  in  the  radius  of 
the  swerve  vector  equal  to  0.015  meters.  This  is  approximately  5  times  the  demonstrated 
range  accuracy  and  results  in  a  Cjva  error  of  about  7%. 

Most  projectiles  fired  supersonically  had  swerve  armi  of  less  than  0.01  meters  and 
expected  coefficient  errors  approach  15%. 

The  multiple  fit  results  are  plotted  in  Figure  3.  The  cubic  term  was  estimated  on 
many  projectiles  and  held  constant.  These  estimates  were  based  on  Reference  1  and  .sev¬ 
eral  sets  of  wind  tunnel  data.  The  accuracy  of  the  zero  yaw  Cvc.  as  represented  by  the 
faired  line  is  approximately  5%  subsonic/transonic  and  10%  supersonic. 


3.  Magnus  Force  Coefficient 

The  Magnus  force  coefficient’s  effect  on  the  motion  of  the  center  of  gravity  is  nearly 
an  order  of  magnitude  smaller  than  the  normal  force  coefficient.  However,  on  those  multi¬ 
ple  fit  groups  where  the  maximum  yaw  observed  exceeded  10  degrees,  this  coefficient  was 
solved  for  and  determined.  The  determined  value  ranged  from  about  -1.0  to  -1.6.  On  all 
other  groups  the  CVpQ  was  estimated  at  -1.0  and  held  constant.  See  Table  3  for  details. 
This  magnitude  of  CVP„  is  similiar  to  the  values  determined  on  other  projectile  test  pro¬ 
grams  with  large  enough  yaw  (Reference  6). 


4.  Pitching  Moment  Coefficient 

The  pitching  moment  coefficient  derivative,  Cma,  is  very  non-linear  with  increasing 
Mach  number  as  shown  in  Figure  3.  This  same  trend  was  first  observed  in  Ref.  1  and 
later  detailed  in  Ref.  2.  The  1987-1989  firings  added  substantially  to  the  definition  of  Cma 
above  M=0.96.  The  Cma  peaks  near  M=0.90  at  a  value  of  4.85.  The  Cma  then  drops  to  a 
low  value  of  4.45  at  M=0.96  and  rapidly  climbs  to  4.63  near  M=1.0.  The  supersonic  peak 
is  near  M=1.5  with  Cma  equal  to  4.88.  A  downward  trend  then  begins  with  Cma  falling  to 
4.51  near  M=2.3.  Tabulated  Cma  are  given  in  Table  3  for  the  multiple  fits.  The  estimated 
accuracy  of  the  zero  yaw  Cma  is  about  2%. 

The  cubic  pitching  moment  is  very  weak  with  a  value  of  about  -7  subsonically  and 
less  than  -3  supersonically.  Estimated  values  were  utilized  in  many  of  the  multiple  fits. 
The  values  were  obtained  from  Ref  2.  and  are  consistant  with  other  experimental  data. 
The  error  in  the  cubic  pitching  moment  coefficient  is  about  +/-  3. 
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5.  Pitch  Damping  Coefficient 


The  pitch  damping  coefficient,  Cmq,  can  only  be  determined  on  spin  stabilized  projec¬ 
tiles  when  the  Magnus  moment  coefficient  is  well  modeled  with  respect  to  non-linearities. 
Linear  Theory  reductions  produce  large  coefficient  scatter  for  Cmq  when  Magnus  non- 
linearities  are  present.  Linear  Theory  values  for  Cmq  varied  from  4-47  to  -30  below  M=l. 
The  6DOF  multiple  fit  analysis  reduced  this  variance  to  a  range  of  -4  to  -17.  Some  of  the 
remaining  variance  is  due  to  Mach  number  trends,  however,  it  is  estimated  that  tv  end  line 
drawn  in  Figure  4  has  an  estimated  error  of  magnitude  3. 

No  attempt  was  made  to  determine  higher  order  terms  of  pitch  damping  as  simulated 
test  cases  have  shown  minimal  influence  on  the  resultant  motion. 


6.  Spin  Damping  Coefficient 

The  spin  damping  coefficient,  C/p,  was  estimated  based  on  yawsonde  data  contained 
in  Reference  1.  This  estimate  is  plotted  in  Figure  4. 


7.  Magnus  Moment  Coefficient 

The  Magnus  moment  coefficient,  Cnp  of  the  M483A1  is  very  non-linear  with  angle  of 
attack  and  Mach  number.  When  the  maximum  yaw  angle  in  the  data  set  to  be  analyzed 
is  greater  than  5  degrees,  the  use  of  higher  order  terms  must  be  closely  reviewed  for 
adequacy.  Reference  2  discusses  the  use  of  an  estimated  Magnus  table  ( Cnpa  vs.  a).  The 
interpolated  derivative  is  then  form  factored  by  using  only  linear  and  cubic  terms  to  evolve 
the  appropriate  non-linear  behavior  with  angle  of  attack. 

For  the  multiple  fit  analysis  of  test  data  below  M=l.l,  a  table  of  Cnp  vs  a  vs  Mach 
number  was  estimated  and  input  in  a  data  file.  As  described  in  the  Data  Analysis  Section 
above,  the  values  of  C„p  were  interpolated  and  then  form  factored  into  a  contour  which 
minimized  the  probable  error  of  fit.  This  process  was  iterative  in  that  the  input  table  was 
adjusted  several  times  to  ensure  realistic  modeling  over  the  entire  range  of  angle  of  attack 
and  Mach  number. 

The  final  coefficients  are  given  in  Table  4.  This  is  the  actual  data  utilized  for  the  final 
multiple  fit  reductions. 

Table  3  (Multiple  Fits)  shows  values  for  Cnpa  above  M=l.l.  Below  M=l.l,  the  ap¬ 
proach  described  above  was  utilized  and  the  expansion  terms  are  not  cubic/quintic.  Figures 
5  thru  16  present  plots  of  Cnp  vs  a  and  Cnpa  vs.  a2  at  Mach  numbers  from  0.65  to  2.0. 

In  general,  below  M=0.98,  Cnpa  is  very  negative  (about  -5)  near  zero  yaw,  crosses  0.0 
at  about  4  degrees  and  increases  to  a  maximum  of  1.0/1.25  at  6  degrees  yaw.  This  trend  is 
similiar  to  that  reported  in  Ref.  2.  However,  the  goodness  of  the  Cnpa  trend  with  yaw  level 
presented  in  Figures  5  to  16  is  much  better  than  Reference  2.  The  large  quantity  of  data 


/ 
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from  M=  0.85  to  M=i.O  greatly  improved  the  quality  of  and  confidence  in  the  results.  The 
probable  errors  of  fit  were  very  low  (most  less  than  0.1  degrees)  considering  that  4  and  5 
projectiles  with  maximum  yaw  levels  exceeding  12  degrees  were  being  simulataneosly  fit. 


8.  Stability  Considerations 

The  gyroscopic  and  dynamic  stability  of  the  M4S3A1  has  been  computed  for  muzzle 
exit  conditions  at  ambient  conditions.  The  equations  given  below  are  discussed  in  Reference 
4. 


Gyroscopic  Stability  Factor 


c  =  _ UlsL _ 

J9  —  npIyCma<PV1 


Ballistic  Factor 


Fast  -  Slow  Vector  Magnitude  (Initial  Conditions) 

I<f  - 

Ks  - 

Mean  Squared  Yaw 

P  =  K\  +  A’J, 

Fast  Mode  Damping  Factor 


Slow  Mode  Damping  Factor 


The  gyroscopic  stability  factor  of  the  M483A1  was  computed  for  muzzle  exit  conditions 
(1/20  twist)  and  both  a  standard  (15C)  and  cold  (-40C)  day.  These  values  are  presented 


in  Table  5. 


The  fast  and  slow  mode  damping  factors  in  the  above  equations  were  set  equal  to 
zero,  the  equations  re-arranged,  and  solved  for  the  two  limit  values  of  Cnpa •  The  slopes 
computed  have  been  plotted  in  Figures  5  thru  16. 

The  point  at  which  Cnp  is  crossed  by  the  negative  slope  is  the  magnitude  of  the  slow 
mode  (arm)  limit  cycle.  The  projectile  will  cone  at  this  yaw  angle  (3.5  degrees  at  M=  0.87 
-  Figure  8)  at  steady  state  conditions. 

If  the  Cnp  is  crossed  by  the  positive  slope,  the  projectile  has  a  fast  mode  dynamic 
instabilty.  Using  Figure  8  (Mach  0.87)  as  an  example,  a  slight  dynamic  instability  caused 
by  a  positive  C„p  exists  from  about  5  thru  9  degrees  yaw.  Above  9  degrees,  the  C„p  falls 
below  the  boundary  and  the  projectile  is  stable  in  both  fast  and  slow  modes. 

The  same  answers  can  be  achieved  using  the  Cnpa  versus  mean  squared  yaw  presen¬ 
tations.  This  trend  was  noted  in  Reference  2. 

The  consequences  of  this  highly  non-linear  dynamic  stability  situation  does  not  readily 
lend  itself  to  closed  form  prediction.  A  series  of  6DOF  trajectory  simulations  with  initial 
conditions  (pitch-yaw  rate)  being  varied  for  purposes  of  computing  the  effect  on  range 
precision  should  be  done.  Cold,  ambient,  and  hot  atmospheres  should  be  included  along 
with  the  effect  on  temperature  on  the  muzzle  velocity. 


IV.  Conclusions 

The  aerodynamic  coefficients  of  the  M483A1  projectile  have  been  computed  based  on 
spark-range  tests  conducted  in  1975,1987,1988,  and  1989. 

The  multiple  fits  showed  that  there  were  essentially  no  aerodynamic  differences  be¬ 
tween  the  projectiles  manufactured  in  the  mid  1970’s  and  late  1980’s. 

The  probable  errors  of  the  multiple  fits  were  very  low  and  support  a  high  confidence 
level  in  the  determined  coefficients.  Overall,  the  quality  of  the  Transonic  Range  data  from 
the  1988  tests  is  the  best,  within  the  experience  of  the  author.  This  indicates  that  both 
the  calibration  of  the  facility,  and  the  film  reading  process  have  improved  in  recent  years. 

The  coefficients  presented  in  this  report  should  be  incorporated  into  an  aerodynamic 
data  package  for  the  M483A1.  This  package  should  be  made  available  to  all  government 
and  private  industry  parties  involved  in  the  development/evaluation  of  projectiles  with 
shapes  similiar  to  the  M483A1. 
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Table  2.  Single  Fit  -  Six  Degree  of  Freedom  Results 
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Table  2.  Single  Fit  -  Six  Degree  of  Freedom  Results 


6  DOF  Suuiry  output 
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Table  2.  Single  Fit  -  Six  Degree  of  Freedom  Results 


|  Probable  Error 

Shot  Buaber  Busbar  DBSQ  CZ  CHa  CTpa  Caa  Caq  Capa  Clp  CXH  j  X(a)  Angle (deg) 

ABABM  CX2  CHa 3  CZpa3  Caa3  Caq2  Cnpa3  Cnpa5  IZ/1Z  CaaM  |  Y-X(a)  Koll(deg) 
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Table  4.  Magnus  Moment  Cnp  vs  Angle  of  Attack  vs  Mach  Number 
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Table  5.  Muzzle  Gyroscopic  Stability  Factor 
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Figure  1.  Configuration  for  155mm  projectile,  M483A1 
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6  DOF  Reduction  -  Coefficient  Output 


Mach  Number 

Figure  2.  Axial  Force  Coefficient  vs  Mach  Number 
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6  DOF  Reduction  -  Coefficient  Output 
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Figure  3.  Normal  Force  and  Pitching  Moment  Coefficients  vs  Mach  Number 
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6  DOF  Reduction  -  Coefficient  Output 
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Figure  5.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  0.65 
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Figure  0.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  0.75 
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M483A1  MACH  NO.  0.870 
ALPHA  VS.  CNP 
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Figure  8.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  0.87 
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M483A1  MACH  NO.  0.900 
ALPHA  VS.  CNP 


- 

M483A1  MACH  NO.  0.900 
OBAR-SQ  VS .  CNPA 


Fast  Mode  Boundary 


Slow  Mode  Boundary 


DELBAR  SQ 

Figure  9.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  0.90 
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Figure  10.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  0.93 
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M483A1  MACH  NO.  0.960 
ALPHA  VS.  CNP 
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ALPHA  VS.  CNP 
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Figure  12.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  1.00 
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M483A1  MACH  NO.  1.250 
ALPHA  VS.  CNP- 


DELBAR  SQ 

Figure  14.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  1.25 
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M483A1  MACH  NO.  1.500 
ALPHA  VS,  CNP 


ALPHA 
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OBAR-SQ  VS.  CNPA 


0  25  50  75  100  125  150  175 


DELBAR  SQ 

ure  15.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  1.50 
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M483A1  MACH  NO.  2.000 
DBAR-SQ  VS .  CNPA 
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Figure  10.  Magnus  Moment  Coefficients  vs  Angle  of  Attack  at  Mach  No.  2.00 
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